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Abstract 



The recent discovery of metadislocations in some periodic complex 
metallic alloys and of their 'phason' defects has given a new impetus to 
the study of QC approximant defects. In this paper we emphasize: 1— 
that approximants differ from a QC by a suitable density of 'flips' only, 
one per unit cell in the case of Fibonacci approximants; these flips are 
not topological defects, 2— a flip can split into two 'phason' defects of 
opposite signs, thus the approximant defects can be studied in a first 
step as defects of the parent QC. In a companion paper this analysis of 
QC defects is extended to dislocations; the difference between perfect 
and imperfect dislocations is emphasized. Imperfect dislocations are 
the phason defects alluded to above. 

1 introduction 

My renewed interest in the quasicrystal defect theory starts with the 
discovery by the Julich group of a family of very remarkable line defects 
present in periodic complex metallic alloys, metadislocations, which display 
an extremely small Burgers' vector compared to the unit cell parameter and 
are attended by specific 'phason' defects, see Ref. [1] for a review. These 
investigations have motivated a reappraisal of the question of approximants 
and of their defects, see e.g. [2, 3], a question previously investigated in 
many papers, [4, 5, 6, 7] to cite a few. It is this author contention that 
metadislocations enter the framework of quasicrystal defects, without need 
for new concepts in the defect theory of the solid state, at least in a first step. 
This paper (I) and the following (II) ref. [8], which constitute an extended 
synthesis of previous published articles [5, 9], are devoted to the definition 
of approximants as quasicrystals with specific defects (I) and to new consid- 
erations on perfect b|| and imperfect bj_ dislocations in quasicrystals (II). 

The sites of a rational approximant are the intersections of a set of atomic 
surfaces (i.e. copies of the dj_-dimensional projection onto Ej_ of a cell of the 
(i-dimensional hypercubic lattice £ ) attached to a d||-rational cut E with an 
irrational (i||-space E|| that carries the physical structure; E, E||, Ej_ C £, 
d = dn + dj_. In this way the resulting atomic distances are the same as in 
the i-phase. In [3] rational approximants are analytically described as the 
result of an homogeneous shear of the underlying hyperspace, as in [4] . 
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We have proposed in [5] a different way of generating an approximant, 
which has the advantage of showing directly how the periodic lattice inherits 
the characteristics of the defects of the parent quasicrystal. The unit cell of 
an approximant is a piece of the parent QC with a certain number of flips, in 
fact only one flip in the simplest ones — the Fibonacci approximants — which 
can conveniently be located at the vertices of the unit cell. We recall that 
a flip is a local displacement of an atom, which, in the tiling representation 
of QCs, consists in the shift of a vertex which respects the tile shapes but 
breaks the tile matching rules; see examples hereunder. This is developed 
in Sect. 2. 

In this paper we present a complete view of the nature of rational ap- 
proximants, how they relate to the parent QCs by the presence of 'flips', and 
how those flips can split into 'phason' defects, which are nothing else than 
the 'imperfect' dislocations of the QC, generally loosely attached to the per- 
pendicular component of the Burgers vector bj_. Flips are not topological 
defects, so that approximants can be thought of as elastic instabilities of a 
parent QC; this is indeed the description adopted by Jaric & al. [4]. On 
the other hand flips can split into 'phason' defects of opposite signs, each 
of them being a true topological defect. Therefore if ever this splitting is 
favored, and opposite phason defects repulsive, the approximant might be 
more stable than the parent QC. A detailed study of the 'phason' defects, 
which are true imperfect dislocations, is made in a companion paper (de- 
fects in quasicrystals, revisited II— perfect and imperfect dislocations), in 
relation with the perfect dislocations b|| to which they can be paired. 

2 construction of a rational approximant 

In this paper we employ the classical construction of a QC: the atoms {m} 
are the intersections with an irrational cut E|| C £ of a set of atomic surfaces. 
We consider in turn the cases: d\\ = d± = 1; du =2, d± = 3; du = d± = 3. 

2.1 Fibonacci approximants in one dimension 

A Fibonacci approximant is represented in a 2D hyperspace by a cut E 
whose slope is the ratio of two consecutive Fibonacci numbers f n , fn+l- A 
rational line of slope p = f n / fn+i is the best approximation to the irrational 
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Fi gure 1 Three periods of a Fibonacci approximant {/3 , /4} in the hyperspace d = 2 (sawtooth 
cut Eg). . . .a + /3~, /3 + 7 _ -y + 8~ . . . are irrational segments parallel to En, slope t~ 1 , linked by 
the 'phason' shifts . . . ct~ a + , (5~ f3 + , f~ 7+ . . . perpendicular to En; ... A, B,C . . . are vertices of 
the lattice. AB = BC = . . . = {/4, 73}. The approximant generated by the sequence of irrational 
cuts is the same as the approximant generated by E (dashed rational line), slope fz/ fA = 2/3. 

line En of slope r _1 (r = (v5 + l)/2 is the golden ratio) in the sense 
that there are no vertices of the square lattice in a double triangle like 
B/3 + (I)7~C, if B and C are vertices of this lattice, see Fig. 1; the segment 
BC has slope p, the segment /3 + 7~ has slope r . This result is independent 
of n; besides the area a of the double triangle above is also independent of 
n, a = 5~2 ; see appendix. 

This has important consequences: a 'phason' shift which moves af5 to 
AB does not meet any vertex of the hyper lattice, and thereby there is no flip 
of the atomic sites during this move. This allows for the representation of 
the cut E of a rational Fibonacci approximant by a sawtooth-like sequence 
Eg of segments of slope t , linked by 'phason' shifts a~a + , f3~ f3 + , ■ ■ ■ , see 
Fig. 1 . Thereby as far as we are interested in the construction of a Fibonacci 
approximant along E, Eg is equivalent to E. 

The approximant can be constructed as follows. Attach an atomic sur- 
face (a copy of the projection of a cell on the perpendicular space) to each 
cell center, Fig. 2, upper left corner; those atomic surfaces that intersect 
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Eg determine a sequence of S (for short) segments of length ? and L (for 
long) segments of length t = re;, when projected upon a copy of E|| . In 
the case of a Fibonacci approximant, the period in 'real' 1-space consists in 
fn+2 = fn + fn+l segments, namely f n Ss and / n+1 Ls. 
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Figure 2: The same as Fig. 1, with some atomic surfaces represented. See text. 

However this process, when applied to the cut Eg (or E) of Fig. 1, pro- 
duces an ambiguity, shown Fig. 2. Because of the high symmetry of the 
Eg + E arrangement with respect to the hyperlattice, the atomic surfaces 
attached to the cells C3 and C3" both do not intersect Eg and E (or possibly 
both intersect Eg and E for different atomic surfaces), so that a choice has 
to be made. It is equivalent to slide Eg to the right (then C3 is active) 
or to the left (then C3" is active ) by a sufficient amount. The sequences 
of active atomic surfaces are then either . . . C2C3 C 4 ... or . . . 020^04 . . . , 
i.e. . . . SL. . . or . . . LS. ... In fact, whatever the choice that is made, the 
sequence of S and L is the same, the difference amounting to a shift in the 
chosen origin of the period; it is . . . LSLLS. . . for any sawtooth-cut Eg with 
n = 3, wherever it is in the hyperlattice. 

This ambiguity being straightened out, we now show that each period 
exhibits a flip with respect to a perfect quasi-lattice; the comparison which 
makes sense is between the sequences carried by E|| and those carried by 
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E (or E5). For this purpose let us compare a sequence of active atomic 
surfaces belonging to E: 

n n n+n n n n n+n 

which yields a sequence of S & L segments: 

LSLLSLSL, 

and those of E|| continued along its segment f3j': 

Cc^ r^~^ — 
1^2^3 W^5^6W^8 *^9> 

which yields a sequence: 

LSLLSLLS. 

There is therefore a flip associated to this transformation E|| -4 E5 along 
the segment BC. More generally the 'phason' shifts aA, /3B, • • • , introduce 
one flip per period, comparing with a local quasilattice carried by the seg- 
ments AB, BC, .... One can consider that the two lattice vertices at the 
two extremities of the triangles, each counting for a J vertex, amount to the 
presence of 1 vertex per double triangle. Thereby it is appropriate to say 
that there is one flip, and only one, per period of a Fibonacci approximant; 
notice furthermore that a appears as the natural area that has to be swept 
by the cut line E|| in a 'phasonic' move that produces one flip, and only one. 

The term flip refers to the fact that two neighboring segments along the 
Fibonacci sequence are interchanged. In terms of an atomic description, 
if such neighboring segments are a S and a L segments, the atom separat- 
ing them has been moved by a distance £ — q = t~ 2 £ along the d\\ = 1 
approximant. 

2.2 Penrose tiling approximants, d = 5, du =2, d± — 3 

Entin-Wohlman & al. [5] have shown, using the pentagrid definition of 
a Penrose tiling, that there is also only one flip in a Fibonacci approximant 
of such a tiling, this flip being anywhere in the unit cell of the approximant. 
This unit cell can be chosen as a parallelogram of edges 

bp = /p+iLp + /pSp, 



bq = fq+l'Lq + fqS q , (1) 
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for which p = q ± 1, or a rhombus for which p = q. Here Lj (resp. Sj) are 
the long (resp. short) diagonals of elementary Penrose tiles, which are either 
thin (resp. thick) rhombi of angles a = 36° (resp. a = 72°). 1 We have 

bp • b q 

— - = cos a, 
|b P ||b 9 | 

a = Z{Lp,Lq} (or a = Z{S P , S g } , S p is parallel to L p , S g to L 9 ). 

In 5-space, this unit cell belongs to a d\\ = 2 rational cut E, with periods: 

P = {0, fp+l, fp, -fp, -fp+l}, q = {—f q +l,0,fq+l, fq,— fq}, (2) 

for a thin unit cell, and 

P = {0)/p+l>/p>— fp>— fp+l}, <l = {—fq,—fq+U®,fq+l,fq}: (3) 

for a thick unit cell. 

As in the previous case where du = 1 , one can use instead of this rational 
cut an irrational sawtooth-cut Eg which is sketched Fig. 3 for the case p = q. 
The rational cut E and the irrational cut E|| both belong to the 4 dimensional 
subspace £4 € £5, orthogonal to the five-fold axis {1,1,1,1,1}, of equation 
J2i x t = 7 = 0, 1 = 1, ... ,5; thus all the operations belong to the same class 
of local isomorphism (LI) in the sense of [11]. We recall that other LI classes 
(other values of 7) correspond to Penrose tiling with different (and more 
complex) arrowing rules [12]. The generic case 7 7^ will not be considered 
here. 

The 4D sawtooth-cut is constructed, by generalizing the 2D case, as 
follows. In £4, E and any copy of E|| intersect, generically, only in one 
point. We choose the vertices of the approximant in the E plane to be 
such intersections. Consider the vertex A in Fig. 3; E|| in A provides an 
element of E5 limited by a double triangle S'~ j3'~ {A)5 + /3 + : in A<5 + /3 + , 5 + 
and /3 + are the projections of B and D on the copy of E|| in A — thus B/3 + 
and D<5 + are 'phason' displacements — ; in A5'~ f3'~ , 5'~ and j3'~ are the 
projections of B' and D' on the same copy of E|| — B'/3 /_ and D'S' being 
also 'phason' displacements. The 4D elements of the sawtooth-cut meeting 
in A are the simplices ABD5 + /3 + and AB'D'<5' /?' with a common vertex 
A, which generalize the 2D elements in Fig. 2, the double simplex ( a double 
triangle in 2D) B/3 + (I)7~C, with I in common. 

1 We recall that these tiles adjust side by side according to matching rules represented 
by a fixed arrowing of each type of tile [10], see hereunder Fig. 4. 
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Figure 3: Geometric relations between Eg, formed of irrational double triangular regions based 
on Ei | , and the rational cut E. E, which is 2D, and its rhombi-made tiling, is 'correctly' sketched 
in 2D. Each 2D copy of E|| intersects E in only one point (4D geometry), here at vertices of the 
rhombic tiling, A, B, C, D, . . . , and is sketched as a thick segment limited to the projections of 
<5, /3 neighboring vertices: 5 + , f} + are the projections of B, D on the irrational plane attached to 
A; <5~, f3~ are the projections of B, D on the irrational plane attached to C. 

As above, there is no vertex of the hypercubic lattice inside the 4D 
volumes of the simplices ABD5 + /3 + and AB'D'J' f3' . Each vertex being 
common to four rhombi in the E plane, there is in fact 4x \ = 1 vertex in 
the ABCD rhombus, and by the same extension as in the ID case, 1 vertex 
per rhombus, thus one flip per unit cell in the Fibonacci approximant of a 
Penrose tiling. 

The flips of a Penrose tiling are well-known objects. We revisit them in 
the next section, in relation with the notion of stacking fault, i.e. a Penrose 
rhombus edge that does not obey the matching rules. 

2.3 3D approximants 

The same arguments can be developed for 3D approximants of i-phases, 
employing 6D simplices that form a sawtooth. We do not have constructed 
such geometries, but there is no doubt they exist, and that they yield Fi- 
bonacci approximants with one flip per unit cell. A phason singularity is 
now a face of a Mackay rhombohedron that does not obey the matching 
rules, in relation with the presence of flips. 

3 about phason defects 

A flip is the sum of two singularities, each of them akin to a stacking 
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Figure 4: Flips a — > a' and b — > b' accompanied by two 'phason' singularities (or mismatches) 
of opposite signs. The two singularities can diffuse apart along the row of hexagons by a sequence 
of local flips. Adapted from Fig. 3 in [9]. 

fault bordered by an imperfect dislocation dipole. Let us summarize what 
this statement implies: 

— In fig. 4 are sketched two flips in a Penrose tiling that bring the atoms 
in positions a, b to the positions a', b'. Because a flip can be erased by a 
movement of the atoms bringing them back to their previous positions, the 
two singularities (the marked edges on the figure) that accompany it can be 
thought of as being of opposite signs; it is equivalent to refer to a flip or to 
these two singularities. The marked edges do not satisfy the matching rules 
— here represented by the de Bruijn's arrowing [10] — in a Penrose tiling 
of isomorphism class 7 = 0. These mismatches are the singularities we are 
alluding to, 

— the two mismatches carried by the flip can diffuse apart [13] by a 
succession of flips, independently one from the other. Therefore it makes 
sense to consider a mismatch as a 'phason' singularity per se, since it cannot 
be eliminated by a local movement of the atoms, 

— mismatches are true stacking faults in the sense of the condensed mat- 
ter defect theory of defects. The Penrose tiling case is illustrated fig. 5. 
The shift that characterizes the fault is a displacement b* equal to the flip, 
fig. 5c). Figures 5 a), b) make explicit the construction of the dislocation 
dipole Li, L/2, from which one deduces its Burgers vector b. The shift b* 
is a part of b by virtue of this construction; it is this shift that constitutes 
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a) b) c) 

Figure 5: Matching fault in a Penrose tiling 7 = 0. A phason singularity is akin to an imperfect 
dislocation dipole: a) construction of a dislocation of Burgers vector L1L2 (b = {1,0,0,0,0} in 
the hyperspace), with a zigzaging cut surface along Li«; after matter removal (blue) the two 
zigzag lines Li», L2/3 are glued. The plastic deformation thus introduced (not drawn) is relaxed 
by a dislocation of opposite Burgers vector in L2, implying matter addition (red), b) Final result: 
L1L2 is a stacking fault (arrowing undefined) bordered by an imperfect dislocation dipole: b docs 
not belong to the set of b Burgers'vectors. c) The dotted lines show the perfect configuration in 
contact with the matching fault, before any flip action; b* is the shift of the vertex that defines 
the matching fault; it is a part of b. Adapted in part from Fig. 4 in [9]. 



the imperfect feature of the dislocation. This will be more thoroughly dis- 
cussed in (II). In 3D, matching faults are made of a 2D fault (an imperfect 
dislocation), 

— stacking faults (also called matching faults) are therefore singularities 
per se, and as such can be classified by the methods of the topological theory 
of defects [14, 15], 

4 discussion 

There is some confusion in the literature about the concept of phasons. 
We have to distinguish between flips and topological defects. A flip is not a 
topological defect, since an antiflip (exerted on the same atom or cluster) 
restores the perfect structure. The topological defects we are alluding to are 
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stacking faults, with their usual meaning; because they have been thought 
of in the wake of the notion of mismatches — a concept that is typical of 
QCs — they can also be called matching faults. I have used both terms, 
without making a distinction. Two opposite topological defects can collapse 
to form a flip. 

The fact that one finds experimentally the same b||S in the approximants 
and the parent QC [1] justifies the representation of an approximant in the 
hyperspace by an saw-tooth-like irrational cut rather than by a continuous 
rational cut. 

Also, since b|| measures a distance between atomic species in a Mackay 
cluster, one can infer that the cluster shape and scale in the QC and in the 
approximant are the same. This was anticipated in Entin-Wohlman & al. 
[5]. The question therefore arises whether this can be given a meaning in 
terms of electronic stability rules, applied respectively to the QC and to the 
approximant [16]? In a sense an approximant might be nothing else than a 
modification of a QC at practically constant energy. Notice that one does 
not know of any faceted crystal of an approximant, whereas quasicrystals 
often grow faceted. 
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appendix A 

Calculation of the length /3~/3 + . 

We use the Binet's formula, namely 

fnV5 = [r n -(-l) n T- n ]. 

The angle w = Zf3~j~ (3 + = v — u, where u and v are the angles of the 
directions E and E| | with the abscissa, namely tan u = /n//n+l> tanu = r _1 . 
Thus: 

tanu> = (tan v — tan it) /(l + tan u tan v) reads, using the formula above: 

tanu; = (-l) n+1 r- 2n - 1 . 
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From that expression one gets sin 2 w = r 2n 1 /(/2n+i\/5)- Since ((3 7 ) 2 = 
Un + /n+i) sin2 w = hn+i sin 2 one eventually gets: 

Calculation of the area of the triangle /3~7~/3 + . 

This is a right triangle, and the double of its area is a = (/3~/3 + ) x 
(/3 + 7~), hence: 

/• 1 1 , I tantu| 1 . 

O = J2n+1 sin W COS W = / 2 n+i— - — ^— = -7= « 0.447214, 

1 + tarn w V5 

which is independent of n. The Binet's formula is useful in the course of the 
demonstration of this result. 
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